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Abstract 

Magnetic resonant coupling (MRC) is an efficient method for realizing the near-field wireless power transfer 
(WPT). Although the MRC enabled WPT (MRC-WPT) with a single pair of transmitter and receiver has been 
thoroughly studied in the literature, there is limited work on the general setup with multiple transmitters and/or 
receivers. In this paper, we consider a point-to-multipoint MRC-WPT system with one transmitter delivering wireless 
power to a set of distributed receivers. We aim to introduce new applications of signal processing and optimization 
techniques to the performance characterization and optimization in multiuser WPT via MRC. We first derive closed- 
form expressions for the power drawn from the energy source at the transmitter and that delivered to the load at 
each receiver. We identify a “near-far” fairness issue in multiuser power transmission due to receivers’ distance- 
dependent mutual inductance with the transmitter. To tackle this issue, we propose a centralized charging control 
algorithm to jointly optimize the receivers’ load resistance to minimize the total transmitter power drawn while 
meeting the given power requirement of each individual load. For ease of practical implementation, we also devise 
a distributed algorithm for the receivers to adjust their load resistance independently in an iterative manner. Last, 
we characterize the power region that constitutes all the achievable power-tuples of the loads via controlling their 
adjustable resistance. In particular, we compare the power regions without versus with the time sharing of users’ 
power transmission, where it is shown that time sharing yields a larger power region in general. Extensive simulation 
results are provided to validate our analysis and corroborate our study on the multiuser MRC-WPT system. 
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I. Introduction 

Inductive coupling |E1, 0 is a traditional method to realize the near-held wireless power transfer (WPT) for short- 
range applications in e.g., centimeters. Recently, magnetic resonant coupling (MRC) lffl-0 has drawn significant 
interest for implementing the near-held WPT due to its high power transfer efficiency as well as long operation 
range, say, up to a couple of meters. Furthermore, MRC effectively avoids the power leakage to non-resonant 
externalities and thus ensures safety to the neighboring environment. 

Two different methods are commonly adopted in practice to implement MRC enabled WPT (MRC-WPT). In 
the hrst method 01, 0, resonators, each of which is a tunable RLC circuit, are placed in close proximity of the 
electromagnetic (EM) coils of the energy transmitters and receivers to efficiently transfer power between them. 
Since resonators are designed to resonate at the system’s operating frequency, the total reactive power consumption 
in the system is effectively minimized at resonance and hence high power transfer efficiency is achieved over 
longer distance as compared to conventional inductive coupling. In the second method f6], Q, series and/or shunt 
compensators, each of which is a capacitor of variable capacity, are embedded in the electric circuits of energy 
transmitters and receivers with their natural frequencies set same as the system’s operating frequency to achieve 
resonance. Generally speaking, the second method achieves higher power transfer efficiency over the first method, 
since in the first method resonators incur additional power loss due to their parasitic resistance. However, the electric 
circuits of energy transmitters and receivers need to be accessible in the second method to embed compensators in 
them. 

The MRC-WPT system with a single pair of transmitter and receiver has been extensively studied in the literature, 
with the aims such as maximizing the end-to-end power transfer efficiency or maximizing the power delivered to 
the receiver’s load with a given input power lf8l- ITO . Moreover, systems with two transmitters and a single receiver 
or with a single transmitter and two receivers have been studied in flT2l - lfl6l . while their results cannot be directly 
applied to the systems with more than two transmitters/receivers. Recently, an MRC-WPT system with multiple 
transmitters and one single receiver has been investigated in fT71 to wirelessly charge a cellphone located at 40 
centimeters away, independent of the phone’s orientation. However, the interactions between the energy transmitters 
and receiver were demonstrated only through simulations in ltl7l . There have been other recent works (see e.g. fl8l . 
||l9l ) on optimizing the performance of MRC-WPT systems with multiple transmitters and one single receiver. 

Different from the above works, in this paper we consider a point-to-multipoint MRC-WPT system based on the 
series compensator method aforementioned, as shown in Fig. |T] where one transmitter that is connected to a stable 
energy source supplies wireless power to a set of distributed receivers. Each receiver is connected to an electric load 
via a switch, where the switch connects/disconnects the load to/from the receiver. We aim to apply signal processing 
and optimization techniques to the performance characterization and optimization in multiuser MRC-WPT systems. 
First, by extending the results in lfT2l - lfI71 . we derive closed-form expressions for the power drawn from the energy 
source at the transmitter and that delivered to the load at each receiver, in terms of mutual inductance among the 
transmitter and receivers as well as their circuit parameters, for arbitrary number of receivers. Our obtained results 
reveal a near-far fairness issue in multiuser wireless power transmission, similar to its counterpart phenomenon in 
multiuser wireless communication. Specifically, a receiver that is far from the transmitter and thus has a small mutual 
inductance with the transmitter receives lower power as compared to a receiver that is closer to the transmitter, 
with other circuit parameters given identical. Next, we propose a method to mitigate the near-far issue by jointly 
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Fig. 1. A point-to-multipoint magnetic resonant coupling enabled wireless power transfer system with communication and control. 


designing the load resistance of all receivers to control their received power by exploiting the mutual coupling 
effect in the MRC-WPT system. This is analogous and yet in sharp contrast to the method of adjusting antenna 
weights at the transmitter to control the received power at different receivers in the existing far-held microwave or 
radio frequency (RF) transmission enabled WPT lf20ll . li2Tl . 

In particular, we consider the scenario where a central controller is equipped at the transmitter to coordinate 
the multiuser power charging, by assuming that it has the full knowledge of all receivers, including their circuit 
parameters and power requirements. The central controller jointly designs the adjustable load resistance of all 
receivers to minimize the total power consumed at the transmitter subject to the given minimum power requirement 
of each load. For ease of practical implementation, we also consider the scenario without any central controller 
installed and devise a distributed algorithm for multiuser charging control by adjusting the loads’ resistance at 
their individual receivers in an iterative manner. In our proposed distributed algorithm, each receiver sets its load 
resistance independently based on its local information and a one-bit feedback broadcasted by each of the other 
receivers. The feedback of each receiver indicates whether the received power of its load exceeds the required 
minimum power level or not. It is shown via simulations that the distributed algorithm achieves performance fairly 
close to the optimal solution by the centralized algorithm with a finite number of iterations. 

Last, we characterize the power region for multiuser power transfer which constitutes all the achievable power- 
tuples for the receiver loads via controlling their adjustable resistance in given ranges. Specifically, we introduce 
the time-sharing based multiuser power transfer, where the transmission is divided into orthogonal time slots and 
within each time slot only a selected subset of receivers are scheduled to receive power, while the other receivers 
are disconnected from their loads. This is aimed to more flexibly control the mutual coupling effect between the 
transmitter and receivers in WPT. It is shown that time sharing can enlarge the power region over the case without 
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time sharing in general. It is also shown that time sharing can further mitigate the near-fare issue in multiuser 
WPT by allocating more time to receivers that are more far-away from the transmitter. Furthermore, we extend 
the centralized multiuser charging control algorithm for the case without time sharing to jointly optimize the time 
allocation and load resistance for all the receivers in the case with time sharing, to further reduce the transmitter 
power consumption under the same average power requirement of each load. 

The rest of this paper is organized as follows. Section II introduces the system model. Section III presents 
our analytical results. Section IV presents both the centralized and distributed multiuser power charging control 
algorithms. Section V characterizes and compares the power regions without versus with time sharing. Finally, we 
conclude the paper in Section VI. 

II. System model 

As shown in Fig. |T] we consider an MRC-WPT system with a single transmitter and N > 1 receivers, indexed 
by n, n £ A f = {1,... , A r }. The transmitter and receivers are equipped with EM coils for realizing wireless 
power transfer, while an embedded communication system is assumed to enable information exchange among 
them [] The transmitter is connected to a stable energy source supplying sinusoidal voltage over time given by 
v tx (t) = Rc{(; tx e r/ ' !/ }, with n tx denoting a complex voltage which is assumed to be constant, and w > 0 denoting 
its operating angular frequency. Each receiver n is connected via a switch to a given electric load (e.g., battery 
charger), named load n, with adjustable resistance x n > 0. The switch is used to connect/disconnect each load 
to/from its corresponding receiver. The state of switch at each receiver n is given by s n € {0, 1}, where s n = 1 
and s n = 0 denote the switch is closed and open, respectively. It is also assumed that the transmitter and each 
receiver n are compensated using series capacitors with capacities c tx > 0 and Cn ^ 0, respectively. 

Let z tx (f) = Rc{i tx e r " ,,: }, with complex-valued z tx , denote the steady state current flowing through the transmitter. 
This current produces a time-varying magnetic flux in the transmitter’s EM coil, which passes through the EM coils 
of nearby receivers and induces time-varying currents in them. We denote i n (t) = Rc{i n e J1 ' :l }, with complex-valued 
i n , as the steady state current at receiver n. It is worth pointing out that the magnetic flux is the main medium 
of wireless power transfer considered in this paper, while the electric field is evanescent and thus is ignored (4). 
This is in contrast to the RF based far-held WPT ll20l . lITTl . where the synchronized oscillations of magnetic and 
electric fields radiate energy in the form of EM waves propagating through the air. 

We denote r tx > 0 (r n > 0) and l tx > 0 (l n > 0) as the internal resistance and the self-inductance of the 
EM coil of the transmitter (receiver n), respectively. We also denote the mutual inductance between EM coils 
of the transmitter and each receiver n by a real number h n , with \h n \ < \//,,/ tx , where its actual value depends 
on the physical characteristics of the two EM coils, their locations, alignment (or misalignment) of oriented axes 
with respect to each other, the environment magnetic permeability, etc. For example, the mutual inductance of two 
coaxial circular loops that lie in the parallel planes with separating distance of d meter is shown to be proportional 
to d~ 3 in ||231 . Moreover, since the receivers usually employ smaller EM coils than that of the transmitter due to 
practical size limitation and they are also physically separated, we ignore the mutual inductance between any pair 
of the receivers for simplicity. 

'As an example, the alliance for wireless power (A4WP) specification (23 uses a low energy profile Bluetooth network at the band 
of 2.4GHz for communication and system control, which is aimed to schedule the charging sequence of receivers and also control their 
individual charging power according to the given priorities. 
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The equivalent electric circuit model of the considered MRC-WPT system is also shown in Fig. |T] in which the 
natural angular frequencies of the transmitter and each receiver n can be expressed as u/ n aturai,tx = 1/v^txCtx ar| d 
^natural ,n = 1 / \Jln ( -'rn respectively. We thus set the capacities of compensators’ capacitors as 


Ctx = 


Cn — 


Ztx W 2: 


1 


In. W‘‘ 


Vn € AT, 


( 1 ) 

( 2 ) 


so that the transmitter and all receivers have the same natural angular frequency as the transmitter voltage source’s 
angular frequency w, i.e., uWurai.tx = '^'natural : i = ... = '(/’natural..v = w. Accordingly, we name w as the resonant 
angular frequency. 

In this paper, we assume that the transmitter and all receivers are at fixed positions and the physical characteristics 
of their EM coils are a priori known. As a result, h n ’s, Vn € A f, are modeled as given constants, which are computed 
according to Appendix E. In practical systems with mobile receivers, h n ’s in general change over time and thus 
need to be measured periodically. For example, one method that can be used in practice to estimate the mutual 
inductance between the transmitter and any receiver n, is given as follows. First, by disconnecting the loads at all 
other receivers k / n, under a known input voltage v tx , the transmitter measures the power drawn from its voltage 
source, denoted by p tx , due to load n only. From ([8]), we can show 


h n = ± 



- Ax (r n + x n ) 


(3) 


w 

i.e., the transmitter can obtain the mutual inductance with receiver n by assuming known r n and x n (which can 
be sent to the transmitter via one-time feedback from receiver n). Note that the sign of h n can be determined by 
comparing the known direction of the current flowing in receiver n (via a one-bit feedback from receiver n) with 
that assumed at the transmitter. If the directions are same, then the positive sign is selected for h n ; otherwise, the 
negative sign is set. 

In this paper, we treat the load resistance x n ’s, Mn e M, as design parameters, which can be adjusted in real 
time to control the performance of our considered MRC-WPT system based on the information shared among 
different nodes in the system, via the embedded communication system. Note that an electric load with any fixed 
resistance can be connected via a rectifier in parallel with a boost (or triboost) converter to each receiver to realize 
an adjustable resistance 11241 . Specifically, given the fixed input voltage, the on/off time intervals of the converter 
can be controlled in real time to change the average current flowing into the load, which is equivalent to adjusting 
the load resistance. 


III. Performance Analysis 

In this section, we present new analytical results on the performance of the MRC-WPT system with arbitrary 
number of receivers. A numerical example is also provided to validate our analysis and draw useful insights. Here, 
we assume that all receiver switches are closed, i.e., s n = 1, Vn € A r ; as a result, the transmitter sends wireless 
power to all loads concurrently. 
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A. Analytical Results 


By applying Kirchhoff’s circuit laws to the electric circuit model given in Fig. [TJ we have 


\r tx + j (wl tx - —) ) i tx - jw V h k i k 

V V wcJJ ^ 


- Ax> 


1 


WC n 


r n + x n + j win - in-jwhnitx = 0, Vn € Af. 


(4) 

(5) 


From ([1} and ©, we can set wl tx — 1 /(wc tx ) = 0 and wl n — l/(wc n ) = 0 in © and I©, respectively. This is 
due to the fact that the transmitter and all receivers are designed to resonate at the same angular frequency w. By 
solving the set of 1 inear equations given in © and ©, we can derive i tx and i n ’s as functions of the input voltage 
v tx as follows: 


fix — 


nx + w 2 Efce./V h l ( r k + x k ) 1 


_1 V tXl 


In =3- 


whn (r n + X n ) 


-1 


—n tx , Vn € Af. 


( 6 ) 


(7) 


rtx + w 2 J2kGAf h l i r k + Xk) 

The power drawn from the energy source at the transmitter, i.e., p tx , and that delivered to each load n, denoted by 
p n , are then obtained as 


Ptx = -Re{u tx i*} 

_ 1 I . i2 _ 

Pn — 2 \^n\ — 


2 r tx + w 2 h l ( r k + x k) 1 ’ 

|n tx | 2 w 2 h 2 n x n \r n + x n )~ 2 _ 

2 (r tx + w 2 Efce M h l ( r k + Zfc) -1 ) 


(8) 

(9) 


where i* x denotes the conjugate of / tx . From ©. it follows that the power delivered to each load n increases with 
the mutual inductance between EM coils of its receiver and the transmitter, i.e., h n . This can cause a near-far 
fairness issue since a receiver that is far from the transmitter generally has a small mutual inductance with the 
transmitter; as a result, its received power is lower than that at a receiver that is closer to the transmitter (thus has a 
larger mutual inductance). Furthermore, we define p sum = Pk as the sum-power delivered to all loads, where 

it can be verified from © and ® that p sum < p tx . The sum-power transfer efficiency, denoted by 0 < p < 1, is 
thus expressed as 


= Psum = ^ 2 Efc€A r h l x k(r k +X k ) 2 1Q) 

Ptx r tx + w 2 J2 kGA f h 2 k (r k + x k )~ l 

Remark 3.1: When the receivers are all weakly coupled to the transmitter, e.g., they are sufficiently far away 
from the transmitter, we have h n —>• 0, Vn E Af. In this regime, from ®, it follows that the transmitter power is 
Ptx ~ |r' tx | 2 /(2r tx ), which is a function of the resistance and voltage of the transmitter only. On the other hand, 
from ©. it follows that the power delivered to each load n is p n ~ |u tx | 2 w 2 h^ l x n {r n + x n )~ 2 /(2r 2 x ), which is 
irrespective of the other receivers’ mutual inductance and resistance. The above results can be explained as follows. 
With h n —>• 0, Vn € Af, the power transfered to the receivers is small and thus can be neglected as compare to the 
power loss due to the transmitter’s resistance. As a result, we have p tx ~ r tx |i tx | 2 /2, with i tx = v tx /r tx . It also can 
be verified that with h n —» 0, Vn € Af, the coupling effect among the receivers through the transmitter current i tx 
is negligible. Hence, the power delivered to the load at each receiver is independent of other receivers (similar to 
the far-field RF based WPT [[20]], EH ). 
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Remark 3.2: It can be shown from © that p n , Mn € N, first increases over 0 < w < w, and then decreases 
over w > w, where 


w = 


Ax 


(ID 


V Efce M K i r k + x k ) 1 

The above result can be explained as follows. From ©, it follows that the magnitude of the current flowing in 
each receiver n, i.e., \i n \, strictly increases over 0 < w < w, but strictly decreases over w > w. This yields that 
w = w is the unique maximizer of \i n \ over w > 0. Obviously, p n , which is defined in © as p n = x n \i n \ 2 /2, 
behaves same as \i n \ over w > 0. Although w is assumed to be fixed in this paper, it can also be optimally set 
to maximize the system power transfer efficiency, if this is implementable in practice. Furthermore, from (fTIT) . it 
follows that w depends on the distances between the transmitter and receivers, since h n ’ s in general decrease with 
larger distances. 

Next, we study the effect of changing the load resistance of one particular receiver n, i.e., x n , on the transmitter 
power ptx, its received power p n . the power delivered to each of the other loads m € AT, m / n, i.e., p rn , the 
sum-power delivered to all loads p sum , and the sum-power transfer efficiency p, assuming that all the other loads’ 
resistance is fixed. 

Proposition 3.1: p lK strictly increases over x n > 0. 

Proof: Please see Appendix A. ■ 

This result can be explained as follows. From ©. it is observed that the transmitter current |i tx | strictly increases 
over x n > 0. Since the energy source voltage u tx is fixed, it follows that p tx given in © strictly increases over 

Xn P D* 


Proposition 3.2: p m . Mm f n, strictly increases over x n > 0. However, p n first increases over 0 < x n < x n , 
and then decreases over x n > x n , where 


r n (Ax + <t>n) + w 2 h 2 n 
Ax T 4*n 


( 12 ) 


with f n = w 2 EfeeAAfn} h K r k + Xk) l - 

Proof: Please see Appendix B. ■ 

The above result can be explained as follows. From ©, it follows that for each receiver m, m 7 ^ n, its current 

\i m \ strictly increases over x n > 0. The received power p rn defined in © thus strictly increases over x n > 0. On 

the other hand, it follows from © that for receiver n, its current \i n \ strictly decreases over x n > 0. However, 

from ©, it follows that the decrement in \i n \ 2 is smaller than the increment of x n when 0 < x n < ,x n ; therefore, 

p n increases over x n in this region. The opposite is true when x n > x n . 

Proposition 3.3: If r tx + f n -2p n < 0, p sum strictly increases over x n > 0, where ip n = w 2 JfkeAT\{n} hl x k(xk + 

Xk)~ 2 ', otherwise, p sum first increases over 0 < x n < x n , and then decreases over x n > x n , where 

„ _ r n (r tx + <j> n ) + w 2 h 2 n + 2r n p n ,, ^ 

X n — 1 „ • O-H 

Ax T fn 2 p n 


Proof: Please see Appendix C. ■ 

This result is a consequence of Proposition 13.21 from which it is known that p m ’s, m n, strictly increase over 
x n > 0, while p n first increases over 0 < x n < x n and then decreases over x n > x n . The sum-power p sum can 
thus behave similarly as either p m ’ s (mono tonic ally increasing) or p n (initially increasing and then decreasing) over 
x n > 0 , depending on the system parameters. 











Fig. 2. The considered system setup for numerical examples. 

Proposition 3.4: If tp n ~ 4>n ~ r i\ > 0 , p strictly increases over x n > 0; otherwise, p first increases over 
0 < x n < x n , and then decreases over x n > x, where 

x = -AM (14) 

Pn tlx 

with r n = {tp n -4> n - r tx )(r2(r tx + p n + 4> n ) + r n w 2 h 2 ). 

Proof: Please see Appendix D. ■ 

This result is a consequence of Propositions 13.11 and 13.31 due to the different behaviors of p tx and p sum over 

Tn T 0. 

B. Validation of Analysis 

For the purpose of exposition, we consider a point-to-multipoint MRC-WPT system with N = 3 receivers, as 
shown in Fig. [2] where the transmitter and receivers use circular EM coils (see Fig. |T2] in Appendix [E]), with 
the physical characteristics given in Table U Note that the transmitter and both receivers 2 and 3 lie in the plane 
with z = 0, while receiver 1 lies in the plane with z = 0.91 meter (m). Accordingly, the internal resistance and 

self-inductance of individual EM coils as well as the mutual inductance among them can be derived (see the details 

in Appendix [E]). where the obtained values are given in Table |TTJ In this example, although all receivers use EM 
coils with the same physical characteristics, they arc located in different distances from the transmitter. Specifically, 
receiver 1 is closest to the transmitter and thus has the largest mutual inductance with the transmitter, while receiver 
3 is farthest and has the smallest mutual inductance. We set v u = 20\/2V, and w = 42.6 x 10 6 rad/s (i.e., 6.78MHz), 
as suggested in the A4WP specification 1221 . For this example, we fix X 2 = %3 = 2.5f2. 

First, we plot p tx , p n ’ s, Vn E J\T, and p sum versus the resistance of load 1, i.e., x i, in Fig. [3] It is observed that 
Ptx, P'2, P 3 and p sum all increase over x\ > 0; however, p\ initially increases over 0 < x\ < i;i = 5.3512 and then 












TABLE I 

Physical characteristics of EM coils 
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EM Coil 

Inner radius 

(cm) 

Outer 

radius (cm) 

Average 

radius (cm) 

Number of 

turns 

Material of 

wire 

Resistivity of 
wire (pQ/m) 

Transmitter 

19.9 

20.1 

20 

200 

Copper 

0.0168 

Receiver 1 

4.95 

5.05 

5 

10 

Copper 

0.0168 

Receiver 2 

4.95 

5.05 

5 

10 

Copper 

0.0168 

Receiver 3 

4.95 

5.05 

5 

10 

Copper 

0.0168 


TABLE II 

Electrical characteristics of EM coils 


EM Coil 

Internal resistance 

f'x.yJfn (^) 

Self-inductance l t Jl n 

(mH) 

Mutual inductance h n 

(pH) 

Transmitter 

1.3440 

54.0630 

- 

Receiver 1 

0.0672 

0.0294 

-0.0921 

Receiver 2 

0.0672 

0.0294 

0.0402 

Receiver 3 

0.0672 

0.0294 

0.0245 


declines over x\ > 5.3517. Note that in this example, the condition r tx + <i> n — 2 ip n < 0 holds in Proposition 13.31 
The obtained results are all consistent with our analysis in Section III-A. Besides, we observe that varying x\ not 
only changes p\, but also the power delivered to other loads. For instance, receiver 1 can help receivers 2 and 3, 
which are farther away from the transmitter, to receive more power by increasing its load resistance x\. This is a 
useful mechanism that will be utilized later in this paper to mitigate the near-far issue. 

Second, in Fig. ED we plot the sum-power transfer efficiency p as a function of x±. It is observed that p follows 
a single-peak pattern over x n > 0, i.e., it first increases over 0 < x\ < x\ = 0.95Q, and then smoothly declines 
over x\ > 0.9517. This result can be verified from Proposition 13.41 by considering the fact that in this example, the 
condition ip n — 4> n — Gx < 0 holds. Note that when x \ — » oo, it follows from © and © that i\ -X 0 and p\ — » 0. 
This is equivalent to disconnecting load 1 from receiver 1, i.e., setting .s'i = 0. As a result, the efficiency converges 
when x\ —>• oo, while the converged value depends on the parameters of the transmitter and the other two receivers. 

Third, we set x\ = X 2 = X 3 = 2.517, and plot the power received by the three loads versus w in Fig. [5] It is 
observed that p\, P 2 , and p:> reach their individual peaks all at w = w = 17.97 x 10 6 rad/sec, which is in accordance 
to Remark lT2l It is worth noting that although w and x n ’s can be jointly designed to achieve better performance 
over the case of optimizing x n ’s only with w being fixed, this problem is challenging to solve and thus is left as 
our future work. 


IV. Multiuser Charging Control Optimization 

In this section, we optimize the receivers’ load resistance x n ’s to minimize the transmitter power p tx subject to 
the given load constraints, by assuming that s n = 1, Vra € AT, i.e., all the receives are connected to their loads. 
First, we consider the case with a central controller at the transmitter, which has the full knowledge of all receivers, 





















10 



Fig. 3. Input and output power versus xi. 



Fig. 4. The sum-power transfer efficiency p versus x\. 



Fig. 5. The output power versus w. 

including their circuit parameters as well as their load requirements, to implement centralized charging control. We 
then devise a distributed charging algorithm for the receivers to independently adjust their load resistance iteratively, 
for the ease of practical implementation. Last, we compare the performance of the two algorithms under a practical 
system setup. 
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A. Problem Formulation 


We assume that in practice the resistance of each load n can be adjusted over a given range x n < x n < x n , 
where x n > 0 (x n > x n ) is the lower (upper) limit of the resistance. It is also assumed that the power delivered 
to each load n needs to be higher than a given minimum threshold p > 0 to guarantee its quality of service. 
Next, we formulate the optimization problem (PI) to minimize the transmitter power p tx subject to the given load 
constraints of all receivers as follows. 

12 i 


(PI) : min 

1 ^ 'Tn } n£J\f 


s.t. 


Ax 


r t x + w 2 J2kGAf h l( r k + x k) 


-l 


Ax 


w h 2 n x n (r n + x n ) 


-2 


' Ax + w 2 Efce M h l ( r k + x k ) 1 

Although (PI) is non-convex, we propose a centralized algorithm to solve it optimally in the next subsection. 


> p , Vn E Af. 


(15) 

(16) 


B. Centralized Algorithm 

First, based on (PI), we formulate the maximization problem (P2), where its objective function is the inverse of 
that of (PI) but with the same constraints as (PI). 

2 


(P2) : 


max 


{x n <x n <x n }neAf Ptx| 


s.t. 


\Vtx\ 


Ax + W 2 ^2 h k {r k + x k ) 1 
V k&Af / 

w 2 h 2 l x n (r n + x n )~ 2 


> p , Vn e J\T. 


(17) 

(18) 


Vtx + W 2 h 2 k (r k + x k ) 1 

It can be verified that the optimal solution to (P2) also solves (PI); as a result, we can equivalently solve (P2) to 
derive the optimal solution to (PI). Although (P2) is still non-convex, we can re-formulate it as a convex problem 
by applying change of variables. Specifically, we define a new set of variables as y n = 1 /{r n + x n ), Vn € M. 
Since x n <x n < x n , it follows that y n <y n < y n , where y = l/(r n +x n ) and y n = l/(r n + x n ). Accordingly, 
we rewrite (P2) as (P3). 

2 


(P3) : 


max 


{y_ n <yn<y n }neM \ v tx\ 

..t.M 


Ax + W 2 ^2 h k Vk 


k&M 


w 2 h 2 n (r n y 2 n - y n ) + p_ n r, x + w 2 ^ h\ y k < 0, Vn E JV. 

V fee Jd ) 


(19) 


( 20 ) 


Note that (P3) is a convex optimization problem, with a linear objective function and linear/quadratic inequality 
constraints over y n ’ s. As a result, (P3) can be efficiently solved using the existing software, e.g., CVX ll26l . Let 
(yl,... ,Vn) denote the optimal solution to (P3). The optimal solution to (P2) is thus obtained by a change of 
variable as x* = 1 /y* — r n , Vn € Af. The obtained (.ij\... ,x* N ) also solves (PI). The centralized algorithm to 
solve (PI) is summarized in Table UTT1 denoted as Algorithm 1. Since the feasibility of convex problem (P3) can 
be efficiently checked, in the rest of this paper, we assume that (PI), or equivalently (P3), is feasible without loss 
of generality. 
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TABLE III 

Centralized algorithm for (PI). 

_Algorithm 1_ 

a) For each receiver n, Vn € A f, given x n > 0 and x n > x n , compute y = l/(r n +Xn) and y n = I /(r„ + x n ). Accordingly, formulate the problem 

(P3). 

b) If (P3) is feasible, then save its optimal solution as (y* * ,...,2/jv)- Set a:* = 1/y* — r n , Vn £ A f. Return (a;*,... ,x* N ) as the optimal solution to 
(PI). 

c) If (P3) is infeasible, then it follows that there is no feasible solution to (PI) and thus the algorithm terminates. 


C. Distributed Algorithm 

In this subsection, we present an alternative distributed algorithm for (PI), for the case without a central controller 
installed in the system. In this algorithm, each receiver adjusts its load resistance independently according to its local 
information and a one-bit feedback received from each of the other receivers indicating whether the corresponding 
load constraint is satisfied or not. We denote the feedback from each receiver n which is broadcasted to all other 
receivers as FB n € {0,1}, where FB n = 1 ( FB n = 0) indicates that its load constraint is (not) satisfied. 

In Section Hill we show that the power delivered to each load n, p n , has two properties that can be exploited to 
adjust x n . First, p n strictly increases over x m > 0, Vm / n, which means that other receivers can help boost p n by 
increasing their individual load resistance. Second, p n has a single peak at x n = x n , assuming that all other load 
resistance is fixed. Thus, over 0 < x n < x n , receiver n can increase p n by increasing x n ; similarly, for x n > x n , it 
can increase p n by reducing x n . Although receiver n cannot compute x n from (fl2l) directly due to its incomplete 
information on other receivers, it can test whether 0 < x n < x n , x n = x n , or x n > x n as follows. Let /> n (x+), 
Pn(x n ), and p n (x~) denote the power received by load n when its resistance is set as x n + Ax, x n , and x n — Ax, 
respectively, where Ax > 0 is a small step size. Assuming all the other load resistance is fixed, receiver n can 
make the following decision: 

• If Pn(x+) > p n (x n ) and p n (x~) < Pn{x n ), then 0 < x n < x n ; 

• If Pn(x+) < Pn{x n ) and Pn{x~) < p n (x n ), then x n = x n B 

• If Pn(x+) < p n (x n ) and p n (x~) > Pn(x n ), then x n > x n . 

Now, we present the distributed algorithm in detail. The algorithm is implemented in an iterative manner, say, 

starting from receiver 1, where in each iteration, only one receiver n adjusts its load resistance, while all the other 
receivers just broadcast their individual one-bit feedback FB m , rn ^ n, at the beginning of each iteration. We 
initialize x n = min{ max{(r n r tx + w 2 li 2 n ) /r tx . x n }, x n }, Vn € A f, where (r n r tx + w 2 h^)/r tx is obtained from d 
by setting (!> n = 0, i.e., assuming that all other receivers have their loads disconnected^! This is a reasonable starting 
point, under which the power delivered to each receiver is maximized (see Proposition 13.21 ). Then, as the algorithm 
proceeds, all receivers can gradually adjust their load resistance to help reduce the transmit power while meeting 
the minimum power constraints of their individual loads. Specifically, at each iteration for receiver n, if p n < p , 
then it will adjust x n to increase p„. To find the direction for the update, the receiver needs to check for its current 
x n whether 0 < x„ < x n , x n = x n , or x n > x n holds, using the method aforementioned. On the other hand, 
if p n > p n , receiver n can increase x n to help increase the power delivered to other loads when there exists any 

2 More precisely, in this case we have x„ — Ax < x„ < x n + Ax. 

’This requires a protocol design so that when each new receiver is added in the system, its mutual inductance h n is measured and x n is 
accordingly computed and initially set at the receiver. 





TABLE IV 

Distributed algorithm for (PI). 
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_ Algorithm 2 _ 

a) Initialize itr = 1 and itr max > 1. Each receiver n sets x n = min{ max{(r n r tx + tn 2 /i^)/r tx , £ n }, x n }. 

b) Repeat from receiver n = 1 to n = N: 

• Receiver n collects FBm from all other receivers m ^ n. 

• Receiver n updates its load resistance x n according to Cases 1-5. 

• If itr = itr m ax , then quit the loop and the algorithm terminates. 

• Set itr = itr + 1. 


m n such that FB m = 0 is received; or it can decrease x n to help reduce the transmitter power when FB m = 1, 
Mm n. In summary, we design the following protocol (with five cases) for receiver n to update x n . 

Case 1: If p n < p n and 0 < x n < x n , set x n 4— min{x n , x n + Ax}. 

Case 2: If p n < p n and x n > x n , set x n •<— max{x n , x n — Ax}. 

Case 3: If p n > p , x n ^ x n , and 3m / n, FB m = 0, set x n <— min{x n , x n + Ax}. 

Case 4: If p n > p , x n ^ x n , and FB m = 1, Vm / n, set x n <— max{x n , x n — Ax}. 

Case 5: Otherwise, no update occurs. 

We set a maximum number of iterations, denoted by ifr max , after which the algorithm will terminate. The above 
distributed algorithm is summarized in TableIIVI denoted as Algorithm 2. It is worth noting that due to the simplicity 
of Algorithm 2 as well as its distributed nature, this algorithm may not converge to the optimal solution to (PI) 
in general, or may even fail to converge to a feasible solution to (PI) in certain cases, as will be shown by the 
numerical example presented next. 

D. Performance Comparison 

We consider the same system setup as that in Section IIII-BI We set x n = Iff and x n = lOOfl, Vn € AT. We also 
set P 1 = P 2 = 17.5W, but vary p 3 over 0 < p 3 < 37.95W, where (PI) can be verified to be feasible in this specific 
region. For Algorithm 2, we use Ax = 10 ~ 3 and itr max = 3 x 10 5 , with itr max > — x fc )/Ax, which is 

sufficiently large such that each receiver n can search for its load resistance x n over the whole range of [x n , x n ] 
before the algorithm terminates. 

Fig. [ 6 ] compares the transmitter power p tx obtained by both Algorithms 1 and 2 versus p.,. In this example, 
Algorithm 2 converges to a feasible solution to (PI) only over 0 < p 3 < 33.75W, while it yields an infeasible 
solution to (PI) if > 33.75W. Particularly, with > 33.75W, the load resistance (xi,X 2 ,X 3 ) obtained via 
Algorithm 2 can satisfy the power constraints of loads 1 and 2, but not that of load 3. Moreover, it is observed that 
with 0 < p 3 < 33.75W, Algorithm 2 achieves almost the same minimum p tx as that by Algorithm 1, which solves 
(PI) optimally. Notice that when p 3 > 33.75W, the obtained p lx via Algorithm 2 is lower than that of Algorithm 
1. However, this result is not meaningful as the solution by Algorithm 2 in this case is not feasible. 

Fig. [7] shows the convergence of Algorithm 2 under the above system setup with ]>., = 30W. It is observed that 
although pi > p ] and p-> > p 0 at the first iteration, we have 7)3 < p,,. As a result, both receivers 1 and 2 help 
receiver 3 (which is most far-way from the transmitter) for receiving more power by lowering their received power 
levels via increasing their individual load resistance. It is also observed that this algorithm takes around 0.4 x 10 5 
iterations to converge, since we use Ax = 10 -3 in the algorithm for updating x n ’s, which is a small step size to 
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Fig. 6. Performance comparison between centralized control (Algorithm 1) versus distributed control (Algorithm 2). 



Fig. 7. Convergence performance of Algorithm 2. 


ensure smooth convergence. In practice, larger step size can be used to speed up the algorithm but at the cost of 
certain performance loss. 


V. Power Region Characterization 

In this section, we characterize the achievable power region for the receiver loads without versus with time sharing. 
First, we propose a time-sharing scheme to schedule multiuser power transfer by connecting/disconnecting loads 
to/from their receivers over time. We then propose a centralized algorithm to jointly optimize the time allocation 
and load resistance of receivers for time sharing based power transmission, by extending that for (PI) in Section Hvl 
for the case without time sharing. Last, numerical examples are provided to compare the power-region performance 
of the multiuser MRC-WPT system without versus with time sharing. 

A. Multiuser Power Transfer with Time Sharing 

As shown in Fig. [8j there are in general Q = 2 v — 1 time-sharing configurations for transferring power to N 
receiver loads depending on the state of each receiver’s switch. We index these configurations by q, q e Q = 
{1,..., Q}. Specifically, let S = {(si,..., sn) \ s n € {0,1}, Vn € AT} denote the set consisting of all possible 
states of receiver switches. Without loss of generality, we remove the trivial case that all switches are open, i.e., 
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Configuration 1 

Configuration 2 

... 

Configuration Q 








s' ^ ^ 

t 2 

t q 


Fig. 8. Time-sharing based multiuser power transfer. 


s n = 0, Vn E A f, from S by setting S <— S \ {(0,0)}. As a result, we have the cardinality of S as |<S| = Q. 
For convenience, we assign an one-to-one mapping between the elements in the two sets Q and S, i.e., we assign 
each configuration q e Q to one switch state (s\,..., sn) £ S. By default, we assign configuration q = 1 to 
(si,..., sat) = (1 ,..., 1), i.e., all switches are closed. 

Let r > 0 denote the total time available for power transmission. Let the time allocated for the power transmission 
under configuration q be denoted by r q , with 0 < r q < r. We thus have Yhq&Q T q — r > where the strict inequality 
occurs when the required energy, p n r, n = 1,.... N, at all loads are satisfied by the end of Q-slot transmissions, 
where the voltage source at the transmitter can be switched off for the remaining time (r — Yh q eQ r <j) > 0 to save 
energyQ Over all configurations, the average transmitter power and the average power delivered to each load n can 
be obtained from ([ 8 ]) and respectively, as follows: 


Ptx = 


E 


I Ax | 


1 


Pn = 


qeQ 2r + w2 K ( r k + x k,q) 

|Ax | 2 W 2 h‘ix n ,q{r n + X n ,q )~ 2 


-1 


E 


9 SS„ 


2 r 


2 '<?’ 


( 21 ) 

( 22 ) 


V, x + w 2 EfceW h l ( r k + x k, q ) 1 

where x n , q is the resistance value of load n under configuration q. Furthermore, M q C J\f denotes the subset of 
receivers with their loads connected under configuration q, while Q n C Q denotes the subset of configurations under 
which receiver n has its load connected. Compared to the previous case without time sharing, the time allocation 
Tq s, Mq € Q, can provide extra degrees of freedom for performance optimization. 


B. Power Region Definition 

The power region is defined as the set of power-tuples achievable for all loads with a given transmission time r 
subject to their adjustable resistance values. Specifically, the power region for the case without time sharing (i.e., 
all loads receive power concurrently) is defined as 

A-without-TS — (^J ( P\ ■ ■ ■ • • P .V) ? (23) 

{x n <X n <Xn}neJf 

with p n ’ s, Vn € A/", are given in ©. Similarly, the power region for the case with time sharing is defined as 

^with-TS = 1 J (pi,...,p N ), (24) 

{0<T„<T,}„ eQ , J2 q£a r Q <T 

with pn s, Vn € A f, are given in (I22T) . It is evident that the power region with time sharing is no smaller than that 
without time sharing in general, i.e., 72 W] th„ui-is C 72 wlt h-Ts 5 since by simply setting n = r and r q = 0, Mq 1, we 

have AVvjUumi-TS — ^-with-TS- 

4 Note that disconnecting all receivers from their loads, i.e. setting s n = 0, Vn G A f, cannot achieve this goal. This is due to the fact that 
the ohmic resistance of the transmitter circuit still consumes power as long as the transmitter voltage source is on. 
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C. Centralized Algorithm with Time Sharing: Revised 


In this subsection, we extend the centralized algorithm in Section ITTl] without time sharing to the case with time 
sharing by jointly optimizing the time allocation and load resistance of all receivers to minimize the transmitter 
power subject to the given load (average received power) constraints. Hence, we consider problem (P4) as follows. 


(P4) : min 

{"tin — ' 5 ' n .9 — } « £ Af, q£ Q q(ZQ 


E 


Ax 


1 


Ax + W 2 EkeK h l i r k + X k,q) 


-1 


-•X 

qeQr, 


|Ax | 5 


W 2 h 2 n X n>q {r n + X ntq ) 


-2 


2 r 


2 ' <? — iu 


(25) 


> p , Vn G A/", (26) 


(27) 


V tx + w 2 £ fcGA/; h 2 (r k + a*,,) 1 
r, < r. 

q&Q 

Although (P4) is non-convex, we can apply the technique of alternating optimization to solve it sub-optimally in 
general, as discussed below. Since (PI) is assumed feasible, the feasibility of (P4) is ensured due to the fact that 
77-without-TS — 77 wl th-TS ■ 

Initialize n = r and x n p = x*, Vn € J\f, where (x *,..., x* N ) denotes the optimal solution to (PI) for the case 
without time sharing. Moreover, initialize r q = 0 and x n>q = min{ max{(r n r tx + w 2 h 2 )/r tx , x n }, x n }, Vn € A f, 
Mq ^ 1. At each iteration itr, itr = 1,2,..., we design r (; ’s and x n . q ’s alternatively according to the following 
procedure. First, we solve (P4) over r (/ ’s, 'iq G Q, while the rest of variables are all fixed. The resulting problem is 
a linear programming (LP) which can be efficiently solved using the existing software, e.g. CVX ll26l . We update 
Tq s as the obtained solution. Next, we optimize the load resistance for different configurations sequentially, e.g., 
starting from configuration 1 to Q. For each configuration q, we solve (P4) over x' n ., ; ’s, Vn € M, with the rest of 
variables all being fixed. We thus consider the optimization problem (P4—q) as follows. 


(P4 - q) : min_ 

{■tin — x ". <J —} " e w 


Ax 


1 


s.t. 


2r r tx + w 2 J2keAr q h l(. r k + x k, q ) 1 9 

I Ax I tn h n x n , q (r n A x U) q) 


Ax + w 2 M h K r k + x k , q 


i-i 


(28) 


T 1 - In ~ Pn,-q> Vn £ ^• ( 2 9) 


where p n _ q = E mG Q„\M |a x\ 2 w 2 h 2 x n!m (r n + x n , m ) 2 r m /(2r(r tx + w 2 J2keAT m h l( r k + x k,m) 1 ) 2 )- For each 
receiver n, its load power constraint given in (l26l ) is re-expressed in (l29l ). where all power terms that do not involve 
x n ,q, Vn G A f, are moved to the right-hand side (RHS) of the inequality, denoted by p n ,- q , which is treated as 
constant in (P4—q). From ( 122h . it follows that p n ,-q denotes the average power delivered to load n under all other 
configurations, m € Q n \ {q\. Problem (P4—q) has the same structure as (PI); as a result, we can solve it using 
an algorithm similar to Algorithm 1. We then update x n ^ s, Vn € A f, as the obtained solution to (P4—q). At 
the end of each iteration itr, we compute p^' 1 as the objective value given in (l25k The algorithm stops when 
Aptx = Pa ' ' * — Pa'' 1 E A p holds, with p^ = oo by default and A p > 0 denoting a given stopping threshold. 
The above alternating optimization based algorithm for (P4) is summarized in Table [Vj denoted as Algorithm 
3. Note that the convergence of Algorithm 3 for (P4) is ensured since the objective value of (P4), i.e., p tx , is 
non-increasing over iterations, while the constraints in (l26l) and (1271) are all satisfied at each iteration. 






TABLE V 

Algorithm for (P4). 
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_Algorithm 3 _ 

a) Initialize itr = 1, Ap > 0 and Ap tx = p® = oo. Initialize T\ — T and x n i = a;*, Vn £ A/”. Initialize T q = 0 and 
x n ,q = min{ max{(r„r tx + w 2 h“^)/r^, x n }, i„}, Vn S AC Vq 1. 

b) While Aptx > Ap do: 

• Solve (P4) over r 9 ’s, Vq 6 Q, assuming that the rest of variables are all fixed. Update r 9 ’s as the optimal solution to the resulting problem. 

• For q = 1 to q = Q do : 

— Solve (P4—q) using similar algorithm as Algorithm 1. Update x n 9 ’s, Vn € As as the solution to (P4—q). 

• Compute 1251 and save the obtained value as p^ r \ Accordingly, set Ap tx = p[^ tr '' — p tx tr T 

• Set itr = itr + 1. 

d) Return x n ,q’s and r q ’s as the solution to (P4). 


D. Numerical Example 

We consider the same system setup as that in Section ITTT-R1 We set x n = Iff and x n = 100Q, Vn S AT. For 
Algorithm 3 in the case with time sharing, we set A p = HP 3 . 

Figs. |9ja),|9jb), and|9tc) show power regions 77 w ith 0 u t-TS versus 77 w j t h-TS given by (l23l) and (1241 ). respectively, under 
three different resonant angular frequencies of w = 14.2 x 10 6 rad/sec, w = 42.6 x 10 6 rad/sec, and w = 127.8 x 
10 6 rad/sec, respectively. For the puipose of exposition, we consider only the two-user case for Fig. [9] by assuming 
receiver 3 is disconnected from its load, i.e., S 3 = 0 (see Fig. 2). It is observed that 72 w ith-TS is always larger than 
77 w ithout-TS> as expected. It is also observed that the power region difference becomes less significant as the operating 
frequency decreases. This result is explained as follows. When w is sufficiently small, the power delivered to each 
load n in the case without time sharing, given in ©, can be approximated as p n ~ | v lx \ 2 w 2 h 2 x n (r n + x n )~ 2 /( 2 r 3 ), 
from which it follows that there is no evident coupling effect among receivers. In this regime, given load resistance 
x n , the power received by load n does not depend on whether the other receivers are connected to their loads or 
not. Thus, time sharing is less effective and hence cannot enlarge the power region over that without time sharing. 
Last, note that since Algorithm 3 for (P4) in general obtains a suboptimal solution, 77 wlt h-Ts shown in Fig. [9] is only 
an achievable power region under the time-sharing scenario. 

Next, we consider again the case with all three users in Fig. 2. We also fix w = 42.6 x 10 6 rad/sec and 71 = p 0 = 
5W. Fig. [TO] compares the transmitter power p tli obtained using Algorithms 1 and 3 over with 0 < 71., < 55.9W, 
where Algorithm 3 takes at most 4 iterations to converge. It is observed that Algorithm 3 achieves lower yi tx than 
Algorithm 1 over all values of j>,,. This result is expected due to the fact that time sharing provides extra degrees 
of freedom for multiuser power transmission scheduling. Consequently, the time allocation and load resistance for 
the receivers can be jointly optimized to further reduce the transmitter power as compared to the case without time 
sharing when only load resistance is optimized. 


VI. Conclusion 

In this paper, we have studied a point-to-multipoint WPT system via MRC. We derive closed-form expressions 
for the input and output power in terms of the system parameters for arbitrary number of receivers. Similar to other 
multiuser wireless applications such as those in wireless communication and far-held microwave based WPT, a near- 
far fairness issue is revealed in our considered MRC-WPT system. To tackle this problem, we propose a centralized 
charging control algorithm for jointly optimizing the receivers’ load resistance to minimize the transmitter power 
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Fig. 9. Power regions without versus with time sharing. 



subject to the given load power constraints. For ease of practical implementation, we also propose a distributed 
algorithm for receivers to iteratively adjust their load resistance based on local information and one-bit feedback 
from each of the other receivers. We show by simulation that the distributed algorithm performs very close to 
the centralized algorithm with a finite number of iterations. Last, we characterize the achievable multiuser power 
regions for the loads without and with time sharing and compare them through numerical examples. It is shown that 
time sharing can help further mitigate the near-far issue by enlarging the achievable power region as compared to 
the case without time sharing. As a concluding remark, we would like to point out that MRC-WPT is a promising 
new research area in which many tools from signal processing and optimization can be applied to devise innovative 
solutions, and we hope that this paper will open up an avenue for more future works along this direction. 


Appendix 


A. Proof of Proposition \3.1 1 
From (HI), it follows that 


<9ptx = Kl 2 _ w 2 h 2 [r n + x n ) 2 _ (30) 

dXn 2 ( r tx + w 2 Y.k=i h l ( r fc + x fc) _1 ) 

where it can be easily verified that dp tx /dx n > 0 over x n > 0. This means that p tx strictly increases over x n > 0. 
The proof of Proposition 13.11 is thus completed. 























































































19 


B. Proof of Proposition \3.2\ 

From ([9]), it follows that for m / n. 


dpm |Ax| 2w x m h m h n (r m T x m ) (r n -F x n ) 


(31) 



(ru + w 2 J2^ =1 hl( r k + x k ) 


where it can be easily verified that dp m /dx n > 0 over x n > 0. This means that p m , rn f n, strictly increases over 
x n > 0. Similarly, for m = n, from ©, it follows that 



(32) 


where <f> n = w 2 J2keAT\{n} ^k( r k + x k ) 1 - It can be easily verified that dp n /dx n > 0 over 0 < x n < x n , with 
x n > 0 given in (fl2l) . and dp n /dx n < 0 over x n > x n . The proof of Proposition 13.21 is thus completed. 

C. Proof of Proposition 13.31 


Since dp &um /dx n = , dpk/dx n , we can easily prove Proposition 13.31 using the same argument for the proof 

of Proposition 13.21 The detail is thus omitted for brevity. 


D. Proof of Proposition 13.41 
From ( fTOb . it follows that 



(33) 


Accordingly, it can be verified that when ip n — 4>n — Ax > 0 , p strictly increases over x n > 0. Otherwise, if 
(p n — fn — Ax < 0 , then p increases over 0 < x n < x n due to the fact that x^{p n — 4> n — r tx ) + 2 r n ip n x n + 
r n w 2 h‘j l + r 2 (pn + d> n + r tx ) < 0 and the rest of terms in the right-hand side of (l33l) are all positive; similarly p 
declines over x n > x. The proof of Proposition 13.41 is thus completed. 

E. Impedance Characterization of EM Coils 

As shown in Fig. |TT] we consider two circular EM coils, indexed by i, i A {1. 2 }, in the free space (no external 
electric and/or magnetic fields exist). Without loss of generality, we assume that the center of EM coil 1 is located 
at the origin, i.e., (x = 0, y = 0, z = 0), and its surface normal vector is given by n\ = z. On the other hand, we 
assume that the center of EM coil 2 is located at {x = x', y = y 1 , z = z') and its surface normal vector is given by 
ri ‘2 = n x px + riypy + n z pz, with \Jn x 2 + n y ,2 + ti 2 2 = E As shown in Fig. Q21 we assume that each EM coil i 
consists of bi closely wound turns of round shaped wire, where the inner radius of coil is denoted by (: mneiv , > 0, 
while the outer radius is denoted by e oute r,i > e; nne r,i- 

Accordingly, the average radius of each EM coil i and the radius of the wire used to build this coil are obtained 

as e ave ,i = (e out er,i + e inne r,i)/2 and e wire ,i = (e oute r,i - e mntTti )/2, respectively. Let n and k denote the resistance 
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Fig. 11. A system of two circular EM coils. 




Fig. 12. Circular EM coil. 


and self-inductance of each EM coil i. Given e w j rc , 7 ; <C e ave ,A i.c., the wire is much thinner than the average radius, 
which is practically valid, we thus have Q: 


n = 




8e a 


k = &fe a ve,tM( ln ( — ) - 2 ), 


(34) 

(35) 


-"\vire,z 


where a t is the resistivity of the wire used in EM coil i and // = Itt x 10 ' N/A 2 , which denotes the mag¬ 
netic permeability of the air. Let h denote the mutual inductance between the two EM coils. By assuming 


d = \Jx /2 + y' 2 + z' 2 e ave ,i, e ave , 2 > i-e-> the distance between the two EM coils is much larger than their 
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average radiuses, we have |23|| : 

7r/x6i6 2 e ave ie ave 2 ^ 

h =--- dcos(t/ ) sm(6' j cos(© Jn x 2 

4CL 1 \ 

+ 3cos(6 )/ ) sin(0 / ) sin(^) / )n 3/! 2 

+ (2cos 2 (0') - sin 2 (6» / ))n 2i2 ^ , (36) 

where O' = cos ~ 1 (z'/d) and 4>' = tan ~ 1 (y'/x'). 
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